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. Summary. Facilitated or kinetically constrained spin models (KCSM) are a class 

of interacting particle systems reversible w.r.t. to a simple product measure. Each 
dynamical variable (spin) is re-sampled from its equilibrium distribution only if the 
■ surrounding configuration fulfills a simple local constraint which does not involve the 

chosen variable itself. Such simple models are quite popular in the glass community 
since they display some of the peculiar features of glassy dynamics, in particular they 
can undergo a dynamical arrest reminiscent of the liquid/glass transitiom. Due to the 
fact that the jumps rates of the Markov process can be zero, the whole analysis of the 
long time behavior becomes quite delicate and, until recently, KCSM have escaped a 
rigorous analysis with the notable exception of the East model. In these notes we will 
mainly review several recent mathematical results which, besides being applicable to 
a wide class of KCSM, have contributed to settle some debated questions arising in 
numerical simulations made by physicists. We will also provide some interesting new 
extensions. In particular we will show how to deal with interacting models reversible 
w.r.t. to a high temperature Gibbs measure and we will provide a detailed analysis 
of the so called one spin facilitated model on a general connected graph. 
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1 Introduction and motivations 

Consider the following simple interacting particle system. At each site of the 
lattice Z there is a dynamical variable a^, called in the sequel "spin", taking 
values in {0, 1}. With rate one each spin attempts to change its current value 
by tossing a coin which lands head with probability p £ (0, 1) and setting the 
new value to 1 if head and if tail. However the whole operation is performed 
only if the current value on its right neighbor is 0. Such a model is known 
under the name of the East model [18] and it is easily checked to be reversible 
w.r.t. the product Bernoulli(p) measure. A characteristic feature of the East 
model is that, when q := 1 — p w 0, the relaxation to the reversible measure 
is extremely slow [12]: 

where Treiax is the inverse spectral gap in the spectrum of the (self-adjoint) 
generator C of the process. Notice that if one writes p = r^-^ then Treiax ~ 

e'^^'^ as /3 ^ oo, a behavior that is refered to as a super- Arrhenius law in the 
physics literature. 

The East model is one of the simplest examples of a general class of inter- 
acting particles models which are known in physical literature as facilitated 
or kinetically constrained spin models (KCSM). 

The common feature to all KCSM is that each dynamical variable, one 
for each vertex of a connected graph Q and with values in a finite set S, 
waits an exponential time of mean one and then, if the surrounding current 
configuration satisfies a simple local constraint, is refreshed by sampling a new 
value from S according to some apriori specified measure v. These models have 
been introduced in the physical literature [19, 20] to model the liquid/glass 
transition and more generally the slow "glassy" dynamics which occurs in 
different systems (sec [32, 10] for recent review). In particular, they were 
devised to mimic the fact that the motion of a molecule in a dense liquid can 
be inhibited by the presence of too many surrounding molecules. That explains 
why, in all physical models, S = {0, 1} (empty or occupied site) and the 
constraints specify the maximal number of particles (occupied sites) on certain 
sites around a given one in order to allow creation/destruction on the latter. 
As a consequence, the dynamics becomes increasingly slow as the density of 
particles, p, is increased. Moreover there usually exist blocked configurations, 
namely configurations with all creation/destruction rates identically equal to 
zero. This implies the existence of several invariant measures (see [26] for a 
somewhat detailed discussion of this issue in the context of the North-East 
model), the occurrence of unusually long mixing times compared to standard 
high-temperature stochastic Ising models and may induce the presence of 
ergodicity breaking transitions without any counterpart at the level of the 
reversible measure [17]. 
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Because of the presence of the constraints a mathematical analysis of these 
models have been missing for a long time, with the notable exception of the 
East model [3], until a first recent breakthrough [12, 13]. 

In this work we partly review the results and the techniques of [12] but we 
also extend them in two directions. Firstly we show that the main technique 
can be adapted to deal with a weak interaction among the variables obtained 
by replacing the reversible product measure with a general high-temperature 
Gibbs measure. Secondly, motivated by some unpublished considerations of 
D. Aldous [2], we analyze a special model, the so called FA-lf model, on a 
general connected graph and relate its relaxation time to that of the East 
model. 



2 The models 

2.1 Setting and notation 

The models considered here are defined on a locally finite, bounded degree, 
connected graph Q = {V, E) with vertex set V and edge set E. The associated 
graph distance will be denoted by •) and the degree of a vertex x by d^;. 
The set of neighbors of x, i.e. y such that d{y,x) — 1, will be denoted 
by Afx- For every subset V c V we denote by dV the set of vertices in 
V \V' with one neighbor in V' . In most cases the graph G will either be the 
d-dimensional lattice Z"* or a finite portion of it and in both cases we need 
some additional notation that we fix now. For any vertex a; € Z'* we define 
the *, thw oriented and the *-oriented neighborhood of x as 

Af* = {yeZ'^: y = x + ^ti "i^i, = ±1, and^i 0} 

= {?; e A^r : y ^ x + Y.t^-^ a^ej, ai > 0} 
1C% = {y e Af* : y = x + Y^^^^ ctie,;, q,; = 1, 0} 

where are the basis vactors of Z''. Accordingly, the oriented and *-oriented 
neighborhoods d+A, d'^A of a finite subset yl c Z"^ are defined as d+A := 
{UxeaJCx} \ A, d^A := {Uxga/C*} \ Tl. A rectangle R will be a set of sites of 
the form 

R := [ai,6i] x • • • x [ad,bd.] 
while the collection of finite subsets of Z"^ will be denoted by F. 

2.2 The probability space 

Let (S, v) be a finite probability space with i^(,s) > for any s S. G d S 
will denote a distinguished event in S, often referred to as the set of "good 
states" , and q = i'(G) its probability. 

Given (S, i') we will consider the configuration space f2 = fly = 
whose elements will be denoted by Greek letters (w, ??...). If 5' = (V, E') is 
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a subgraph of Q and to G Qy we will write for its restriction to V . We 
will also say that a vertex x is good for the configuration uj \i lo^ & G. 

On Q equipped with the natural cr-algebra we will consider the product 
measure ji := Ylxev ^x, fx = f- If Q' = (^') -E") is a subgraph of Q we will 
write [ly or /ig/ for the restriction of /i to fly' ■ Finally, for any / G L^(^), we 
will use the shorthand notation to denote its expected value and Var(/) 
for its variance (when it exists). 

2.3 The Markov process 

The general interacting particle models that will be studied here are Glauber 
type Markov processes in f], reversible w.r.t. the measure /i and characterized 
by a finite collection of influence classes {Cx}xgVj where Cx is just a collection 
of subsets of V (often of the neighbors of the vertex x) satisfying the following 
general hypothesis: 

Hpl For all X gV and all A G Cx the vertex x does not belong to A. 
Hp2 r := sup^. sup^^^ ^(^> ^) < +00. 

In turn the influence classes together with the good event G are the key 
ingredients to define the constraints of each model. 

Deflnition 1. Given a vertex x € V and a, configuration u), we will say that 
the constraint at x is satisfied by co if the indicator 




1 if there exists a set AG Cx such that uiy G G for all y G A 
otherwise 



(1) 

is equal to one. 

Remark 1. The two general hypotheses above tell us that in order to check 
whether the constraint is satisfied at a given vertex we do not need to check 
the current state of the vertex itself and we only need to check locally around 
the vertex. This last requirement can actually be weakened and indeed, in 
order to analyze certain spin exchange kinetically constrained models [11], a 
very efficient tool is to consider long range constraints ! 

The process that will be studied in the sequel can then be informally described 
as follows. Each vertex x waits an independent mean one exponential time and 
then, provided that the current configuration co satisfies the constraint at x, 
the value ujx is refreshed with a new value in S sampled from u and the whole 
procedure starts again. 

The generator C of the process can be constructed in a standard way (see 
e.g. [27, 26]) and it is a non-positive self-adjoint operator on L^{f2,fi) with 
domain Dom{C) and Dirichlet form given by 
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^(/) = II A* (cx Var,(/)) , / e Dom{C) 

Here \&Xx{f) = J du{u}x)f'^{i^) — (/ di'{(^x)f{^))'^ denotes the local variance 
with respect to the variable lo^ computed while the other variables are held 
fixed. To the generator £ we can associate the Markov semigroup Pt := e*^ 
with reversible invariant measure fi. 

Notice that the constraints Cxioj) are increasing functions w.r.t the partial 
order in H for which ui < iv' iS ui'^ & G whenever iVx G G. However that does 
not imply in general that the process generated by jC is attractive in the sense 
of Liggett [27]. 

Due to the fact that in general the jump rates are not bounded away 
from zero, irreducibility of the process is not guaranteed and the reversible 
measure /i is usually not the only invariant measure (typically there exist 
initial configurations that are blocked forever). An interesting question when 
Q is infinite is therefore whether ^ is ergodic/mixing for the Markov process 
and whether there exist other ergodic stationary measures. To this purpose 
it is useful to recall the following well known result (see e.g. Theorem 4.13 in 
[27]). 

Theorem 2.1 The following are equivalent, 

(a) limt^oo Ptf = m ^^(m) for all f e i^(/u). 

(b) is a simple eigenvalue for C. 

Clearly (a) implies that limj^oo A* ifPtg) = l^{f)l^{9) for any f,g e L'^in), i.e. 
H is mixing. 

Remark 2. Even if ii is mixing there will exist in general infinitely many sta- 
tionary measures, i.e. probability measures jl satisfying fiPt = jj, for all t>0. 
As an example, assume not identically equal to one and take an arbitrary 
probability measure p, such that /i({S' \ G}^) = 1. An interesting problem is 
therefore to classify all the stationary ergodic measures fl of {Pt}t>o, where 
ergodicity means that Ptf = f {ft a.e.) for alH > implies that / is constant 
{fl a.e.). As we will see later, when Q = Z'^ and for a specific choice of the 
constraint known as the North-East model, a rather detailed answer is now 
available [26]. 

When Q is finite connected subgraph of an infinite graph Q^o = (Kxj, -^oo), 
the ergodicity issue of the resulting continuous time Markov chain can be 
attacked in two ways. 

The first one is to analyze the chain restricted to a suitably defined ergodic 
component. Although such an approach is feasible and natural in some cases 
(see section 6 for an example), the whole analysis becomes quite cumbersome. 

Another possibility, which has several technical advantages over the first 
one, is to unblock certain special vertices of G by relaxing their constraints 
and restore irreducibility of the c;hain. A natural way to do that is to imagine 
to extend the configuration oj, apriori defined only in V, to the vertices in 
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\ V and to keep it there frozen and equal to some reference configuration 
r that will be referred to as the boundary condition. If enough vertices in 

\ V are good for r, then enough vertices of G will become unblocked and 
the whole chain ergodic. 

More precisely we can define the finite volume constraints with boundary 
condition r as 

c^,y(w) :=Ca:(u;-T) (2) 

where are the constraints for Goo defined in (1) and u ■ t G f2 denotes 
the configuration equal to uj inside V and equal to r in 14o \ V. Notice that, 
for any x GV, the rate cj_y(a;) (2) depends on t only through the indicators 
{1t^^g}zgb, where B is the boundary set B := {Voo \ V)r\{[JzevCz)- Therefore, 
instead of fixing r, it is enough to choose a subset A4 C B, called the good 
boundary set, and define 

c^vi^):=cly{u;) (3) 

where r is any configuration satisfying Tz G G for all 2; e and ^ G for 
z G B\A4 . We will say that a choice of A4 is minimal if the corresponding chain 
in G with the rates (3) is irreducible and it is non-irreducible for any other 
choice A4' C A4. The choice A4 = B will be called maximal. For convenience 
we will write r™**^ ('^™™) for the corresponding generators. 

Rem,ark 3. Without any other specification for the influence classes of the 
model it may very well be the case that there exists no boundary conditions for 
which the chain is irreducible and/or their existence may depend on the choice 
of the finite subgraph G- However, as we will see later, for all the interesting 
models discussed in the literature all these issues will have a rather simple 
solution. 

We will now describe some of the basic models and solve the problem of 
boundary conditions for each one of them. 

2.4 0-1 Kinetically constrained spin models 

In most models considered in the physical literature the finite probability space 
(5*, v) is a simple {0, 1} Bernoulli space and the good set G is conventionally 
chosen as the empty (vacant) state {0}. Any model with these features will 
be called in the sequel a "0-1 KCSM" (kinetically constrained spin model). 
Although in most cases the underlying graph C/ is a regular lattice like Z'', 
whenever is possible we will try to work in full generality. 

Given a 0-1 KCSM, the parameter q = v{Q) can be varied in [0, 1] while 
keeping fixed the basic structure of the model {i.e. the notion of the good 
set and the constraints c^) and it is natural to define a critical value qc as 



qc = infjg G [0, 1] : is a simple eigenvalue of £} 
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As we will prove below coincides with the bootstrap percolation threshold 
qiip of the model defined as follows [34] ^. For any 77 e 1? define the bootstrap 

map T : Q n &s 

{Tr])^ = if either ri^ = or c^{r]) = l. (4) 

Denote by /i^"^ the probability measme on J? obtained by iterating n-times 
the above mapping starting from jj,. As n 00 jj,'^^^ converges to a limiting 
measure fi^°°^ [34] and it is natural to define the critical value q^p as 

qbp = mi{q G [0, 1] : = -^o} 

where So is the probability measure assigning unit mass to the constant con- 
figuration identically equal to zero. In other words qbp is the infimum of the 

values q such that, with probability one, the graph Q can be entirely emptied. 
Using the fact that the c^j's are increasing function of rj it is easy to check 
that = So for any q > qbp ■ 

Proposition 2.2 ([12]) qc = qbp and for any q> qcO is a simple eigenvalue 
for C. 

Rem,ark 4- In [12] the proposition has been proved in the special case Q = TJ^ 
but actually the same arguments apply to any bounded degree connected 
graph. 

Having defined the bootstrap percolation it is natural to divide the 0-1 KCSM 
into two distinct classes. 

Definition 2. We will say that a 0-1 KCSM is non cooperative if there exists 
a finite set B CV such that any configuration rj which is empty in all the sites 
of B reaches the empty configuration (all O's) under iteration of the bootstrap 
mapping. Otherwise the model will be called cooperative. 

Remark 5. Notice that for a non-cooperative model the critical value qc is 
obviously zero since with /^-probability one a configuration will contain the 
required finite set B of zeros. 

We will now illustrate some of the most studied models. 

[1] Prederickson- Andersen (FA-jf) facilitated models [19, 20]. In the 

facilitated models the constraint at x requires that at least j < dx neighbors 
are vacant. More formally 

C^ = {AcK-- \A\ > j} 

When j = 1 the model is non-cooperative for any connected graph G and 
ergodicity of the Markov chain is clearly guaranteed by the presence of at 

' In most of the bootstrap percolation literature the role of the O's and the I's is 
inverted 
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least one unblocked vertex. When j > 1 ergodicity on a general graph is more 

delicate and wc restrict ourselves to finite rectangles i? in Z''. In that case and 
for the most constrained cooperative case j = d among the irreducible ones, 
irreducibility is guaranteed if we assume a boundary configuration identically 
empty on Quite remarkably, using results from bootstrap percolation [34] 
combined with proposition 2.2, when Q = and 2 < j < d the ergodicity 
threshold qc always vanishes. 

[2] Spiral model [9, 8] This model is defined on with the following choice 
for the influence classes 

= {NE^ U SE^; SE,, U SW,^\ SW,^ U NW,^] NW^ U NE,^} 

where NE^ = (x + 62,0; + ei + 62), SE^ ~ (a; + ei, a; + ei — 62), SWx = 
{x — 62, a; — 62 — ei) and NWx = (x — ei; x — ei + 62). In other words the 
vertex x can flip iff either its North-East {NE^) or its South- West {SW^) 
neighbours (or both of them) are empty and either its North- West {NWx) or 
its South-East [SEx) neighbours (or both of them) are empty too. The model 
is clearly cooperative and in [8] it has been proven that its critical point qc 
coincides with 1 — Pg, where p° is the critical treshold for oriented percolation. 
The interest of this model lies on the fact that its bootstrap percolation is 
expected to display a peculiar mixed discontinuous/critical character which 
makes it relevant as a model for the liquid glass and more general jamming 
transitions [9, 8]. 

[3] Oriented models. Oriented models are similar to the facilitated models 
but the neighbors of a given vertex x that must be vacant in order for x 
to become free to flip, are chosen according to some orientation of the graph. 
Instead of trying to describe a very general setting we present three important 
examples. 

Example 1. The first and best known example is the so called East model [18]. 
Here Q = 'L and for every x e Z the influence class Cx consists of the vertex 
X + 1. In other words any vertex can flip iff its right neighbor is empty. The 
minimal boundary conditions in a finite interval which ensure irreducibility 
of the chain are of course empty right boundary, i.e. the rightmost vertex 
is always unconstrained. The model is clearly cooperative but qc = Q since in 
order to empty Z it is enough to start from a configuration for which any site 
X has some empty vertex to its right. One could easily generalize the model 
to the case when is a rooted tree (see section 6). In that case any vertex 
different from the root can be updated iff its ancestor is empty. The root itself 
is unconstrained. 

Example 2. The second example is the North-East model in 1? [25] . Here one 
chooses Cx as the North and East neighbor of x. The model is clearly coop- 
erative and its critical point qc coincides with \ — p", where p" is the critical 
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threshold for oriented percolation in 1? [34] . For such a model much more can 
be said about the stationary ergodic measures of the Markov semigroup P^. 

Theorem 2.3 ([26]) If q < Qc the trivial measure Si that assigns unit mass 
to the configuration identically equal to 1 is the only translation invariant, 
ergodic, stationary measure for the system. If q>qc the reversible measure n 
is the unique, non trivial, ergodic, translation invariant, stationary measure. 

Example 3. The third model was suggested in [3] and it is defined on a rooted 
(finite or infinite) binary tree T. Here a vertex x can flip iff its two children 
are vacant. If the tree is finite then ergodicity requires that all the leaves of 
T are unconstrained. It is easy to check that the critical threshold satisfies 
qc = 1/2, the site percolation threshold on the binary tree. 



3 Quantities of interest and related problems 

Back to the general model we now define two main quantities that are of 
mathematical and physical interest. 

The first one is the spectral gap of the generator C, defined as 

gap(£):= inf -Afi- (5) 
° ^ ' /#const Var(/) ^ ' 

A positive spectral gap implies that the reversible measure is mixing for the 
semigroup Pt with exponentially decaying correlations: 

Var {Ptf) < e-2tsap(£) Var(/), V / e L\ij). 

Remark 6. In the sequel the time scale T^^i '■= gap^^ which is naturally fixed 
by the spectral gap will be refered to as the relaxation time of the process. 

For a 0-1 KCSM, two natural questions arise. 

1. Define the new critical point q'^ := mf{q G [0, 1] : gap(£) > 0}. Obviously 
9c — 9c- Is it the case that equality holds ? 

2. If q[. = qc what is the behaviour of gap(£) as g J, ? 

As we will sec later for most of the relevant models it is possible to answer in 
rather detailed way to both questions. 

The second quantity of interest is the so called persistence function (see 
e.g. [23, 36]) defined by 

F{t) := J diiiv) P(aJ(s) =Vo,ys<t) (6) 

where {cr]}s>o denotes the process started from the configuration rj. In some 
sense the persistence function, a more accessible quantity to numerical sim- 
ulation than the spectral gap, provides a measure of the "mobility" of the 
system. Here the main questions are: 
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1. What is the behavior of F{t) for large time scales ? 

2. For a 0-1 KCSM is it the case that F{t) decays exponentially fast as t — > oo 
for any q > q'^ ^ 

3. If the answer to the previous question is positive, is the decay rate related 
to the spectral gap in a simple way or the decay rate of F{t) requires a 
deeper knowledge of the spectral density of £ ? 

4. Is it possible to exhibit examples of 0-1 KCSM in which the persistence 
function shows a crossover between a stretched and a pure exponential 
decay ? 

Unfortunately the above questions are still mostly unanswered except for the 
first two. 

3.1 Some useful observations to bound the spectral gap 

It is important to observe the following kind of monotonicity that can be 
exploited in order to bound the spectral gap of one model with the spectral 

gap of another one. 

Definition 3. Suppose that we are given two influence classes Co and C'q, de- 
note by Cx(w) and c'^{u)) the corresponding rates and by £ and CJ the associated 
generators on If, for all w & O and all x & V, c^w) < Cx{ijj), we say 

that C is dominated by C . 

Remark 7. The term domination here has the same meaning it has in the 
context of bootstrap percolation. It means that the KCSM associated to C is 
more constrained than the one associated to C 

Clearly, if jC is dominated by C , 'D'{f) < ©(Z) and therefore gap(£') < 
gap(£). 

Example 4- Assume that the graph G has n vertices and contains a Hamilton 
path r = {xi,X2, ■ ■ ■ ,x„}, i.e. d{xi+i,Xi) = 1 for all 1 < i < n — 1 and 
Xi ^ Xj for all i j. Consider the FA- If model on Q with one special vertex, 
e.g. Xn, unconstrained (cx^ = 1). Then, if we replace Q hy F equipped with its 
natural graph structure and we denote by C and C the respective generators, 
we get that gap(£) > gap(£'). Clearly C describes the FA-lf model on the 
finite interval [1, . . . , n] C Z with the last vertex free to flip. This in turn is 
dominated by ^East, the generator of the East model on [1, . . . ,n], which is 
known to have a positive [3, 12] spectral gap uniformly in n. Therefore the 
latter result holds also for gap(£') and gap(£) . 

Example 5. Along the lines of the previous example wo could lower bound the 
spectral gap of the FA-2f model in Z'', d> 2, with that in Z^, by restricting 
the sets A e Co to e.g. the (ei, e2)-plane. 
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For a last and more detailed example of the comparison technique we refer 
the reader to section 6. 

Although the comparison technique can be quite effective in proving posi- 
tivity of the spectral gap, one should keep in mind that, in general, it provides 
quite poor bounds, particularly in the limiting case q I qc- 

The second observation we make consists in relating gap(£) when the un- 
derlying graph is infinite to its finite graph analogue. Fix r G V and let 
Gn,r C be the connected ball centered at r of radius n. Suppose that 
inf„ gap(£g^'^) > 0. It is then easy to conclude that gap(£) > 0. 

Indeed, following Liggett Ch.4 [27], for any / G Dom{£.) with Var(/) > 
pick /„ £ L^i^^, fJ:) depending only on finitely many spins so that fn^f and 
£/„ ^ £/ in L2. Then Var(/„) ^ Var(/) and P(/„) V{f). But since /„ 
depends on finitely many spins 

Var(/„)=Vara„,^(/„) and V{U)=Vg^JU} 

provided that m is a large enough square (depending on /„). Therefore 



Var(/) 

and gap(£) > inf„ gap(£g„^) > 0. 



> inf gap(£g„^) > 0. 



4 Main results for 0-1 KCSM on regular lattices 

In this section we state some of the main results for a general 0-1 KCSM on 
Z'' which have been obtained in [12]. 

Fix an integer length scale £ larger than the range of the constraints and let 
Z'*(^) = iZ'''. Consider a partition of Z*^ into disjoint rectangles := Aq + z, 
z e Z<^(£), where Ao = {xgZ'^ : < Xi < £ - 1, i = 1, .., d}. 

Definition 4. Given e S (0, 1) we say that Gi C {0, 1}"^" is a e-good set of 
configurations on scale £ if the following two conditions are satisfied: 

(a) n{Gi) > 1 - e. 

(b) For any collection {^^^^}xeic^ of spin configurations such that ^^^^ e Gg 

for all X € ICq , the following holds. For any ^ G ^2 which coincides with 
in UxeiC^A^x, there exists a sequence of legal moves inside LixeK^Ag^ 
(i.e. single spin moves compatible with the constraints) which transforms 
^ into a new configuration t Q such that the Markov chain in Aq with 

boundary conditions t is ergodic. 

Remark 8. In general the transformed configuration t will be identically equal 
to zero on O'^Aq. It is also clear that assumption (b) has been made having in 

mind models like the East, the FA-jf or the N-E which, modulo rotations, arc 
dominated by a model with influence class Cx entirely contained in the sector 
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{y : y = X + X^iLi cti^ij cti > 0}. If this is not the case one should instead 
use a non rectangular geometry for the tiles of the partition of Z'', adapted 
to the choice of the influence classes. For example for the Spiral Model the 
basic tile at lenght scale ^ is a quadrangular region TZq with one side parallel 

to ei and two sides parallel to ei + 62, TZo Ll{So + {i — + 62) with 

So := {x G Z'^ : < Xi < £ — 1, X2 = 0}. In this case condition (b) should also 
be modified by substituting everywhere d'^Ao with 9^ylo := ei,ei — 62, — £2- 

With the above notation the first main result of [12] can be formulated as 
follows. 

Theorem 4.1 There exists a universal constant cq S (0, 1) such that, if there 
exists £ and a e^-good set on scale I, then inf^igF gap(£™'*^) > 0. In par- 
ticular gap(£) > 0. 

In several examples, e.g. the FA-jf models, the natural candidate for the event 
Ge is the event that the tile Aq is "internally spanned", a notion borrowed 
from bootstrap percolation [1, 34, 14, 24, 15]: 

Definition 5. We say that a finite set F G If^ is internally spanned hy a 
configuration r/Gfiif, starting from the configuration rj'" equal to one outside 
r and equal to r] inside F, there exists a sequence of legal moves inside F 
which connects rj^ to the configuration identically equal to zero inside F and 

identically equal to one outside F. 

Of course whether or not the set Aq is internally spanned for 77 depends only 
on the restriction of 77 to tIq. One of the major results in bootstrap percolation 
problems has been the exact evaluation of the /i-probability that the box Aq 
is internally spanned as a function of the length scale £ and the parameter q 
[24, 34, 14, 15, 1]. For non-cooperative models it is obvious that for any q> 
such probability tends very rapidly (exponentially fast) to one as ^ ^ 00, 
since the existence of at least one completely empty finite set B + x C Aq (see 
definition 2), allows to empty all Aq. For some cooperative systems like e.g. 
the FA-2f in Z^, it has been shown that for any q > such probability tends 
very rapidly (exponentially fast) to one as f ^ 00 and that it abruptly jumps 
from being very small to being close to one as £ crosses a critical scale ^0(9)- In 
most cases the critical length £c{q) diverges very rapidly as q J, 0. Therefore, 
for such models and £ > £c{q), one could safely take as the collection of 
configurations 77 such that Aq is internally spanned for 77. We now formalize 
what we just said. 

Corollary 4.2 Assume that lim^^oo m(^o is internally spanned ) = 1 and 
that the Markov chain in Aq with zero boundary conditions on UxeiC^Ai^ is 
ergodic. Then gap(£) > 0. 

We stress that for some models a notion of good event which differs from 

requiring internal spanning is needed. This is the case for the N-E and Spiral 
models, as can be immediately seen by noticing that at any length scale it 
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is possible to construct small clusters of particles in proper corners of the 
tiles that can never be erased by internal moves. The choice of the proper 
e-good set of confugurations for N-E has already been discussed in [12]. For 
the Spiral Model the definition which naturally arises from the results in [9] 
is the following. Let TZq be the region obtained from TZq by subtracting two 
proper quadrangular regions at the bottom left and top right corners, namely 
TZq :=TZo\ {TZbi U TZir) where TZm (TZtr) have the same shape of TZq shrinked 
at length scale £/4: and have the bottom left (top right) corner which coincides 
with the one of TZq- The e-good set of configurations on scale £, Ge, includes 
all configurations t] such that there exists a sequence of legal moves inside 
TZq which connects t]^" (the configuration which has all ones outside TZq and 
equals r] inside) to a configuration identically equal to zero inside TZq. Lemma 
4.7 and Proposition 4.9 of [9] prove, respectively, property (a) and (b) of 
Definition 4 (with substituted with ei, ei — 62, —62, see remark 8) when 

the density is below the critical density of oriented percolation. Thus, using 
this definition for the good event and Theorem 4.1 we conclude that 

Theorem 4.3 gap(£gpj^g;) > at any p < p°. 

The second main result concerns the long time behavior of the persistence 
function F{t) defined in (6). 

Theorem 4.4 Assume that gap(£) > 0. Then F{t) < e-«sapt + g-Ps^p*. 

Remark 9. The above theorems disprove some conjectures which appeared 
in the physics literature [21, 23, 5, 6], based on numerical simulations and 
approximate analytical treatments, on the existence of a second critical point 
q'^ > Qc at which the spectral gap vanishes and/or below which F{t) would 
decay in a stretched exponential form ~ exp(— f/r)^ with /? < 1. 

Theorem 4.4 also indicates that one can obtain upper bounds on the spec- 
tral gap by proving lower bounds on the persistence function. Concretely a 
lower bound on the persistence function can be obtained by restricting the 
/z-average to those initial configurations r] for which the origin is blocked with 
high probability for all times s < t. Unfortunately in most models such a 
strategy leads to lower bound on F{t) which are usually quite far from the 
above upper bound and it is an interesting open problem to find an exact 
asymptotic as t — > 00 of F{t). 

Finally we observe that for the North-East model on at the critical 
value q = Qc the spectral gap vanishes and the persistence function satisfies 

dtF{yft) = 00 (see Theorem 6.17 and Corollary 6.18 in [12]). 

4.1 Some ideas of the strategy for proving theorems 4.1, 4.4 

The main idea behind the proof of theorem 4.1 goes as follows. First of all one 

covers the lattice with non overlapping cubic blocks {Aix}^^^d and, on the 
rescaled lattice Z''(£) := iZ'^, one considers the new model with single spin 



14 N. Cancrini, F. Martinelli, C. Roberto^', and C. Toninellit' 

space 5 = {0, 1}^'', good event G := Gi, single site measure the restriction of 
H to S and ronormalizcd constraints {c^^"}xez<i(«) which are a strengthening 
of the North-East ones namely 



Such a model is referred to in [12] as the *- general model. By assumption 
the probability of G can be made arbitrarily close to one by taking £ large 
enough and therefore, by the so called Bisection-Constrained approach which 
is detailed in the next section for the case when is a high temperature 
Gibbs measure, the spectral gap of the *-general model is positive. Next one 
observes that assumption (b) of the theorem is there exactly to allow one to 
reconstruct any legal move of the *-general model, i.e. a full update of an 
entire block of spins, by means of a finite (depending only on i) sequence of 
legal moves for the original 0-1 KCMS. It is then an easy step, using standard 
path techniques for comparing two different Markov chains (see e.g. [33]), 
to go from the Poincare inequality for the *-general model to the Poincare 
inequality for the original model. 

The proof of (a slightly less precise version of) Theorem 4.4 given in [12] 
is based on the Feynman-Kac formula and standard large deviation consider- 
ations. However it is possible to provide a simpler and more precise argument 
as follows. One first observe that F{t) = Fi{t) + Fo{t) where 



and similarly for Fo(t). Consider now Fi{t), the case of Fo{t) being similar, 
and define TA{r}) as the hitting time of the set ^ := {77 : % = 0} starting 
from the configuration 77. Then (see e.g. Theorem 2 in [4]) 



where denotes the probability over the process started from the equilibrium 
distribution ji and Xa is given by the variational formula for the Dirichelt 
problem 



Notice that for any / as above Var(/) > n{A) = q. Therefore > 9 gap and 
the proof is complete. 

4.2 Asymptotics of the spectral gap near the ergodicity threshold. 

An important question, particularly in connection with numerical simulations 
or non-rigorous approaches, is the behavior near the ergodicity threshold qc 
of the spectral gap for each specific model. Here is a set of results proven in 



d;^{ri) = 1 iff ?7y e G for all y G tC*. 







[12]. 
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Eeist Model. 



limlog(l/gap)/(log(l/g))2 = (2 log 2) 

g->0 



-1 



(8) 



FA-lf. For any d > 1, there exists a constant C = C{d) such that for any 

q G (0, 1), the spectral gap on Z'^ satisfies: 



as g I 0, where exp"* ^ denote the {d— l)*'^-iterate of the exponential function 
and Ai = 7rVl8. 

The proof of the lower bounds is a rather delicate combination of the renormal- 
ization scheme described above together with paths techniques as described 
in [33]. The upper bounds arc proved instead cither by a careful choice of a 
test function in the variational characterization of the spectral gap or by a 
lower bound on the persistence function F{t) combined with the upper bound 
given in Theorem 4.4. 

Remark 10. Again some of the above findings disprove previous claims for the 
East model [36] and for the FA-lf model in d = 2, 3 [6]. The result for the 
East model actually came out as a surprise. In [36] the model was considered 
"essentially" solved and the result for the spectral gap was gap « q^°&2{i) as 
g I to be compared to the correct scaling g'°S2(9)/2, \^ [3] ^j^g above solution 
was proved to be a lower bound and an upper bound of the form g'°S2(9)/2 
was rigorously established but considered poor because off by a power 1/2 
from the supposedly correct behavior. 

The scaling indicated in [36] is based in part on the following consideration. 
Fix g <C 1 and consider the East model on the interval Aq := [0, . . . , l/g] with 
the last site free to flip {i.e. zero boundary conditions). Notice that 1/g is the 
average distance between the zeros. Start from the configuration identically 
equal to one and let T be the (random) time at which the origin is able to flip. 
Energy barriers consideration (see [2, 3, 16]) suggest that E(T) should scale 
as g''°S2(9) and that is what was assumed in [36]. However it is not difficult to 
prove that the scaling of E(r) is bounded above by (ggap)~^. Indeed we can 
write for any t> 




exp<'-i(c/g2)] '<gap(£)< [exp'^-i(^l-^)] ' 

exp(-c/g«) <gap(£) < exp(-^^l^) 



rf> 3 



d = 2 



(9) 



exp(-cggap(£yijt) > F{t) > fx{al\ ones)P(r >t)> e"¥(T > t) 
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where F{t) is the finite volume persistence function. Integrating over t 
and using the monoticity of the gap (sec [12, Lemma 2.11]) give E(T) < 
e^c(ggap(£yi^)))~^ < e^c(q'gap(£))~^. This, in view of Theorem 4.1, is in- 
compatible with the assumed scaling q^°S2iQ)_ 

Moreover one can obtain a lower bound on E(T) as follows. Let A be such 
that P(r > A) = e-i then clearly P(r > t) < e-L*/^J and E(T) > e'^X. We 
can always couple in the natural way two copies of the process, one started 
from all ones and the other from any other configuration r], and conclude that 

P(the two copies have not coupled at time t) < P(T >t)< e^"-^*. 

Standard arguments give immediately that gap~^ < A i.e. E{T) > e^^ gap~^. 
In conclusion 

e-i[gap(£Aj]"' < ]E(T) < e^qgapiLAj]'' 
5 Extension to interacting models 

In this section we show how to extend the results on the positivity of the 
spectral gap for 0-1 KCSM on a regular lattice Z"^ to the case in which a weak 
interaction is present among the spins. We begin by defining what we mean 

by an interaction. 

Definition 6. A finite range interaction ^ is a collection # := {<?yi}yieF where 
i) (1>A ■ i^A^ for every A G F; 

ii) = z/diam(yl) > r for some finite r = r(<?) called the range of the 

interaction: 

iii) ||<?|| = sup^gzd XIabx 11*^/11100 < oo; 

We will say that (p G BM,r if r{<P) < r and ||^|| < M. 

Given an interaction G B^^m and yl G F, we define the energy in yl of a spin 
configuration cr G 1? by 

Ha{'j)= ^^(^) 

For IT G Oa and r G f2A<^ we also let H]^{a) :— ifyi((T • r) where a ■ r denotes 
the configuration equal to a inside A and to r outside it. Finally, for any 
yl G F and r G 17^= , we define the finite volume Gibbs measure on with 
boundary conditions r and apriori single spin measure v by the formula 

"^A xeA 

where Z^'^ is a normalization constant. 
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The key property of Gibbs measures is that, for any V C A and any ^ in 

A\V , the conditional Gibbs measure in A with boundary conditions r given ^ 
coincides with the Gibbs measure in V with boundary condition r^c • ^. More 
formally 

^iT{■\<^v^ = = ^^V''■H■) 

Clearly averages w.r.t. 1 (yv- = £.) arc function of ^ and, whenever con- 

fusion does not arise, we will systematically drop ^ from our notation. 

As it is well known (see e.g. [35]), for any r < oo there exists Mq > 
such that for any < M < Mq the following holds. For any € Br^M there 
exists a unique probability measure /i* on J7, called the unique Gibbs measure 
associated to the interaction ^ with apriori measure v, such that, for any r, 

lim u.' = u 

where the limit is to be understood as a weak limit. Moreover the limit is 
reached "exponentially fast" in the strongest possible sense. Namely, for any 
d C yl e F and any two boundary conditions r, r', 



max 



Z^^-^ <i^p,(T,T')|e-'"4d.^^(-.-')j (10) 



Ma (<^d) 



where Dr{T,T') = {y : < d{y,A) < r, Ty ^ r'y} and the constants m,K 
depend only on M,r,d. Moreover m "f +oo as M [ 0. When d = 1 the 
threshold Mq can be taken equal to +oo. In all what follows we will always 
assume that ^ G Br,M for some apriori given r, M and that M < Mq. 

Remark 11. In general the constant Mq does not coincide with any "critical 
point" for the model. It is only a sort of "high temperature threshold" (see 
[29, 30, 31] for more details about this issue). 

Having described the notion of the unique Gibbs measure corresponding to 
we can define the generator of a 0-1 KCSM with interaction ^ and 
constraints Cx given by (1), as the unique self-adjoint operator on L?{Q,^^) 
with quadratic form 

X 

where now the local variance Var*(/) is computed with the conditional Gibbs 
measure given all the spins outside x. The construction of the generator in 
a finite volume A with boundary conditions r is exactly the same as in the 
non-interacting case and we skip it. 



5.1 Spectral gap for a weakly interacting North- East model 



Instead of trying to prove a very gc;neral result on the spectral gap of a weakly 
interacting KCSM, we will explain how to deal with the interaction in the 
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concrete case of the North-East model introduced in section 2.4. Moreover, in 

order not to obscure the discussion with rcnormahzation or block construc- 
tions, we will make the unnecessary assumption that the basic parameter q of 
the reference measure v is very close to one. 

Theorem 5.1 Let {cx}xe7? be those of the North-East model. There exists 
Qo € (0, 1) and for any r < oo there exists Mi such that, for any M < 
min(Mo, Mi) and q> qo, 

inf gap(£*)>0 

Remark 12. As wc will see in the proof of the theorem, the restriction on 
strength of the interaction comes from two different requirements. The first 
one is that the finite volume Gibbs measure has the very strong mixing prop- 
erty uniformly in the boundary conditions given in (10). That, as we pointed 
out previously, is guaranteed as long as M < Mq. The second one requires 
that the zeros, which certainly percolate in a robust way w.r.t. the unper- 
turbed measure v because of the assumption q « 1, continue to do so even 
when wc switch on the interaction. It is worthwhile to observe that for the 
one dimensional East model, the first requirement is satisfied for any M < oo 
and that the second one is simply not necessary. Therefore for the East model 
the above theorem should be reformulated as follows. 

Theorem 5.2 Let {cxjxez he those of the East model. For any finite pair 
{r,M) 

inf gap(£*)>0 

Proof (of Theorem 5.1). We will follow the pattern of the proof for the non 
interacting case given in [12] and we will establish the stronger result 

sup7(yl) < +0O, where ^{A) := ( inf inf gap(£^'^)) (11) 

provided that g > go is large and M is taken sufficiently small. Above Maxyi 
denotes the set of configurations in Oa-^ which are identically equal to zero on 
d^A. In what follows in order to simplify the notation we will not write the 
dependence on the boundary condition of the transition rates. 

As in [12] the first step consists in proving a certain monotonicity property 
of 7(^). 

Lemma 5.3 For any V C Ae¥, 

< 7(F) < 7(/l) < 00 

Proof (Proof of the Lemma). Fix G Br,M and, for any ^ G Maxy, define the 
new interaction as follows: 
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if A n F V 

Ea': A'nV=A ^A' {<JA ■ U'\a) if A C F 



Notice that, by construction, 



r(#«) < r(^) and sup ^ ||^'i||oo < ll^*!! 



oo 



so that e Br,M- Next observe that the Gibbs measure on A with interaction 
is simply the product measure 

^^A {(^a) := A^y'^((Ty) (g) J^/i\v(o-yi\y) on = 

Thus, for any / € L"^ {f2v , IJ^'^) and r € MaxA, we can write (Var*'^ = 



where, in the last inequality, we used the fact that, for any x V and any 
u) e Qa, Cx,/i(w) < Cx,v{i^) because ^ e Maxy, together with 



Thanks to Lemma 5.3 we need to prove (11) only when A runs through all 
possible rectangles. For this purpose our main ingredient will be the bisection 
technique of [28] which, in its essence, consists in proving a suitable recursion 
relation between spectral gap on scale 2L with that on scale i, combined 
with the novel idea of considering an accelerated block dynamics which is 
itself constrained. Such an approach is referred to in [12] as the Bisection- 
Constrained or B-C approach. 

In order to present it we first need to recall some simple facts from two 
dimensional percolation. 

A path is a collection {xo.xi, . . . ,Xn} of distinct points in I? such that 
d{xi,Xi+i) = 1 for all i. A *-pa,th is a collection {.to,.xi, . . . ,a;„} of distinct 
points in I? such that Xi+i € N*. for all i. Given a rectangle A and a direction 
ej, we will say that a path {xq, . . . , a;„} traverses A in the i*'*-direction if 
{a;o, . . . , Xn} C A and xq, a;„ lay on the two opposite sides of A orthogonal to 
ej. 

Definition 7. Given a rectangle A and a configuration lj G a path 

{xo, . . . ,Xn} is called a top-bottom crossing ('left-right crossing^ if it tra- 
verses A in the vertical (horizontal) direction and = for all i = 0, . . . ,n. 
The rightmost (lower-most) such crossings (see [22] page 317) will be denoted 
by IIu 



Var <ji,x) 



Var*'«(/) = Var* < j{A)V^^ < 7{A)Vp^{f) 



Var* I {cTy}yjix) = Var*'^(/ 1 {ay}yjix). 



□ 
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Remark 13. Given a rectangle A and a path F traversing A in e.g. the vertical 

direction, let Ap consists of all the sites in A which are in F or to the right of 
it. Then, as remarked in [22], the event {w : il^^ = F} depends only on the 
variables Wx with x € Ap. 

We are now ready to start the actual proof of the theorem. At the beginning 
the method requires a simple geometric result (see [7]) which we now describe. 

Let Ik ■■= (3/2)''/^ and let Ffe be the set of all rectangles yl C which, 
modulo translations and permutations of the coordinates, are contained in 
[0, ^fe+i] X [0, lk+2\ - The main property of Ffe is that each rectangle in Ffc\Ffe_i 
can be obtained as a "slightly overlapping union" of two rectangles in Ffc_i. 

Lemma 5.4 For all k € for all A G Ffe\Ffe_i there exists a finite sequence 



{A^^\ A2^}^^^ in Ffe_i, where Sk := L^fe^^J' ^'^c/i that, letting dk ■= ^Vh — 2, 

(i) A = Af^uA^\ 
(u) diA\A['\A\4'^)>Sk, 



The B-C approach then establishes a simple recursive inequality between the 
quantity "fk '■= sup^gp^ li^) on scale k and the same quantity on scale k — 1 
as follows. 

Fix A e Ffe\Ffe_i and write it as A = A1UA2 with Ai,A2 G Ffc_i satisfying 
the properties described in Lemma 5.4 above. Without loss of generality we 
can assume that all the horizontal faces of Ai and of A2 lay on the horizontal 
faces of A except for the face orthogonal to the first direction ei and that, 
along that direction, Ai comes before A2. Set d = /li fl /I2 and write, for 
definiteness, d = [ai, 61] x [02, 62]. Lemma 5.4 implies that the width of d in 
the first direction, 61 — ai, is at least Sk- Set also 



and let dri = {bi} x [02, 62] be the right face of I along the first direction. 

Definition 8. Given a configuration u G O we will say that co is I -good iff 
there exists a top-bottom crossing of I. 

Given r € Maxyi, we run the following constrained "block dynamics" on 
(in what follows, for simplicity, we suppress the index i) with boundary 
conditions r and blocks Bi := Ai \ I, B2 := A2. The block B2 waits a mean 
one exponential random time and then the current configuration inside it is 
refreshed with a new one sampled from the Gibbs measure of the block given 
the previous configuration outside it (and r outside A). The block Bi does the 
same but now the configuration is refreshed only if the current configuration 
w in B is /-good (see Figure 5.1). 

The generator of the block dynamics applied to / can be written as 




7 = [ai + (61 - ai)/2, bi] x [02,62] 
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-Bi B2 ~ A2 

Fig. 1. The rectangle A divided into two blocks Bi and B2. The grey region is the 
strip I with a top-bottom crossing. 

/:biock/ - - /) + /4f (/) - / (12) 

and the associated Dirichlet form is 

-DtiUf) = l^T (ci Varl^ (/) + Var|, (/)) 
where Ci{oj) is just the indicator of the event that lu is /-good. 
Remark 14-- The reader should keep in mind that the e.g. the notation 

Var|^(/)) 

stands for Var*'^^(/) and that one can imagine the sum re- 

stricted to those configurations outside Bi that coicide with r outside A since 
otherwise their probability /^^'''(^) is zero. 

In order to study the mixing property of the chain we need the following two 
lemmas. 

Lemma 5.5 ([28]) Fix (r, Af) with M < Mq. Then, for any <P e Br,M, 
sup l/x^;'^ (g) - ^^'^''{g)\ < Afe ||.g||oo V.g : Hb- R 

r' 

sup\l4i2 (9) - ^J'A'^i9)\ < Afc ||.g||oo V.g : Hb-^'R (13) 

where Xk '■= Xr/fc+ie"™**/^ and the constants K,m are given in (10). 

Lemma 5.6 There exists qo € (0, 1) and for any r < 00 there exists Mi such 
that, for any M < min(Mo, Mi) and q > qo, 

Ek ■= max max /x^'^(aj is not I -good ) < 6"**=. 
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Proof. It follows immediately from standard percolation arguments together 
with 

sup &\xpif,'JJa^ = !)<(!- q)e'^'^ 

□ 

Wc can now state the main consequence of Lemma 5.5, 5.6. 

Proposition 5.7 There exists qo G (0, 1) and for any r < oo there exists Mi 
such that, for any M < min(Mo, Mi) and q > qo, 

7^?ik:= sup sup (gap(£^;,k))"' < (l - (14) 

for all k so large that the r.h.s. of (14) is smaller than 2. 

Proof (Proof of the proposition). Fix r, M and as prescribed, let r e Maxyi 

and, in order to simplify the notation, drop all the superscripts ^, t. Let 
/ : f2j\ H^- M be a mean zero function, the eigenvalue equation associated to 
the generator (12) is 

Cl(MB,(/)-/)+MB,(/)-/ = A/ (15) 

By construction A > — 2. 

Assume that A > — 1 + \/Afe since otherwise there is nothing to be proved. 
By applying to both sides of (15) and using (13) we obtain 



(1 + A)msJ = /XBi(mS2(/)) ||MBi(/)I|oo< VAfc||MB.(/)||oo (16) 

If we rewrite (15) as 

/ = -f— TT — l^B^if) + T-TTT — /^Si(/) 

1 + A + ci 1 + A + ci 

and apply fiB2 to both sides, by using (16) together with the assumption 
A > -1 + -v/Afe, we get 

||MB2(/)I|oo < ||AiB2(/)l|oo||A*B2( 



1 + A + Cl ' 



+ Afe|| ||oo||MB,(/)||oo 

1 + A + Cl 

< ||MB2(/)I|oO (|IMB2(y^^^)IIcX> + (17) 

which is possible only if 

\M—^—)\\oo>i-Vx'k 

1 + A + Cl 

i.e. 



A < -l + 8v'2Afc + efe 
and the proof is complete. □ 
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By writing down the standard Poincare inequality for the block auxiliary 
chain, we get that for any / 

Var*'^(/) < jI^U Mr(ci Var|^(/) + Var|,(/)) (18) 

The second term in the r.h.s. of (18), using the definition of jk and the fact 
that B2 = A2 G Ffc_i is bounded from above by 

M!'^(var|^(/)) < 7fc-i ^ I^T {c.,b, Varf (/)) (19) 
xeB2 

Notice that, by construction, for all x G B2 and all u), Cx.BjC^) = Cx,a{'^)- 
Therefore the term X]s6B2 I^a'^ {^x,B2 Varf (/)) is nothing but the contribution 
carried by the set B2 to the full Dirichlct form 2?*'^(/). 

Next we examine the more complicate term /i^^'^ ^ci Var*^ (7)^ For any 
Lo such that there exists a rightmost crossing 77^^ in / denote by the set 
of all sites in A which are to the left of U^. Since Var^^ (/) depends only on 
loa\Bi and, for any top-bottom crossing F of /, l^jj^^p} does not depend on 
the variables w's to the left of F, we can write 

(ci Var|, (/)) = /.^'^ (l^g (Var|^ (/))) (20) 

The convexity of the variance implies that 

Ml(Var|^(/))<Varl(/) 

where it is understood that the r.h.s. depends on the variables in il^j and to 

the right of it. The key observation at this stage, which explains the role and 
the need of the event {3 il^^ in /}, is the following. For any w such that ilt^ 
exists the variance Var*^ (/) is computed with boundary conditions (r outside 
A and ujji\/[^ ) which belong to Maxyi^ . Therefore we can bound it from above 
using the Poincare inequality by 

Varl (/) < 7(^a,)I?l (/) < 7(Si U /)I?1 (/) 

where we used Lemma 5.3 together with the observation that Ai_j C Bi U / = 
Ai. In conclusion 

Mr(l{3ii.in/}Ml(Var|^(/))) 

<7(^l)/x5'^(l{3n.in7}I'l(/)) 

< 7(^i)/x^'"(l{3 77. in 7} E ^^-^^ Varf (/)) 

<7(^i)M$'^(E^-.^Varf(/)) 

xeAi 
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because, by construction, for every w such that there exists in I 

Cx,yl„(w) = Ca;^^(w) ^X&A^. (21) 

If we finally plug (5.1) into the r.h.s. of (20) and recall that Ax G J^k-i^ 
we obtain 

Mr(ciVar|^(/)) <7fe-i/x*'^(E ^-.^Var*(/)) (22) 

In conclusion we have shown that 

Var5'^(/) < 7Sck7.-i(2^*'^(/) +E^^(^-.^Var.(/))) (2^) 

xed 

1 /3 

Averaging over the = [/^ J possible choices of the sets , A2 gives 

Var^(/) < 7Sck7fc-i (1 + f (24) 

which implies that 

k 

7fe < (1 + f )7^Sck7fe-i < 7feo n(l + |)7^Lk (25) 

where ko is the smallest integer such that 7yock < 2. If we now recall the 
expression (14) for 7b'iock together with Lemma 5.5 and 5.6, we immediately 
conclude that the product Hjlfeo 7biock(-'- + 7") bounded. □ 



6 One spin facilitated model on a general graph 

In this section we prove our second set of new results by examining the one spin 
facilitated model (FA-lf in short) on a general connected graph Q = {V,E). 
Our motivation comes from some unpublished speculation by D. Aldous [2] 
that, in this general setting, the FA-lf may serve as an algorithm for informa- 
tion storage in dynamic graphs. 

We begin by discussing the finite setting. Let r be one of the vertices and 
T be a rooted spanning tree of Q with root r. On f2 = {0, 1}^ consider the 
FA-lf constraints: 

Cx,g{ijj) = 1 if Wj, = for some neighbor y oi x , . 

otherwise 

and let Cx,g = Cx,g ii x r and 0^,5 = 1. Let £ be the corresponding Markov 
generator and notice that associated Markov chain is ergodic since the vertex 
r is unconstrained. For shortness we will refer in the sequel to C as the {Q, r, 
FA-lf) model. Our first result reads as follows. 
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Theorem 6.1 



gap(a,r,FA-lf) > gap(Z,East) 



Proof. By monotonicity Cj:,g{uj) > Cx.t{'-^) and therefore gap(C/. r, FA-lf) > 
gap(T, r, FA-lf). We can push the monotonicity argument a bit further and 
consider the following (T, r, East) model: 



if either x 
otherwise 



r ov u>y = where y is the ancestor (in T) of x 



(27) 

Clearly Cx,t{co) > Cx,r{ijj) and therefore gap(5, r, FA-lf) > gap(T, r, East). 
We will now proceed to show that 



gap(T, r, East) > gap(Z, East) 



(28) 



If all the vertices of T have degree 2 with the exception of the root and the 
leaves, i.e. if T c Z, then (28) follows from [12, Lemma 2.11]. Thus let us 
assume that there exists x & T with d^; > 3 and let us order the vertices 
of T by first assigning some arbitrary order to all vertices belonging to any 
given layer (= same distance from the root) and then declaring a; < y iff 
either d{x,r) < d{y,r) or d(x,r) = d{y,r) and x comes before y in the order 
assigned to their layer. Let v be equal to the root if d^ > 2 or equal to the first 
descendant of r with degree d^ > 3 otherwise and let Fy = {r,vi, . . . ,Vk,v} be 
the path in T leading from r to v. Let a be a child of v and let 7^ = {Va, Ea) be 
the subtree of T rooted in a. Finally we denote by A and B the two subgraphs 
of T: A := A U T„, B:=T\ %. (see Fig 6). 




' B 



Fig. 2. The subtrees A and B 



Lemma 6.2 

gap(T, r, East) > min (gap(A, r, East), gap(i?, r. East)) (29) 
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By recursively applying the above result to A and B separately, we imme- 
diately reduce ourselves to the case of a tree T' C Z and the proof of the 
theorem is complete. 
□ 

Proof (of Lemma 6.2). In i^(i7,/i) consider the set Tig of functions / that do 
not depend on Ux, x € 7^. Because of the choice of the constraints Cx,t(w), 
Hb is an invariant subspace for the generator of the (T, r, East) model and 

'""U =gap(B,r,East) (30) 

Let us now consider the orthogonal subspace Wg. Any zero mean element 
/ € satisfies i-tTa,{f ) = and therefore we can write 

Var(/) = /x(VarrJ/)) < /x(VarA(/)) 

< gap(^, r, East)"^ ^ /x {c^^a Vara,(/)) 

xeA 

< gap(A, r, East)"^2?(/) 

where the first inequality follows from convexity of the variance and the second 
one is nothing but the Poincare inequality for the East model in A. The proof 
of the Lemma follows at once from (30), (6). □ 

Theorem 6.1 has two consequences that will be the content of the following 
Theorems. The first one deals with the case of an infinite graph. The second 
one deals with the FA-lf model on general graph Q without the special un- 
blocked vertex r but with the Markov chain restricted to a suitable ergodic 
component. 

Theorem 6.3 Let be an infinite connected graph of bounded degree and 
let C be the generator of the FA-lf model on Goo with constraints {cx,g^, x € 
Voo}, i-e. no apriori unblocked vertex. Then 

gap(5oo, FA-lf) > gap(Z,East) 

Proof. The proof combines Theorem 6.1 together with the finite subgraph 
approximation described in section 3. □ 

Theorem 6.4 Let Q be as in Theorem 6.1 and let he the FA-lf generator 
with constraints {cx,g}xi^v on the restricted configuration space i7+ := {77 S 
J7 : Ylxi^vi^ ~ ^a;) — -*-} ^'U'ipP^d with the reversible measure /ti+ := /i(- 1 J7+). 
Then 

gap(£+) > igap(Z,East)M(l2+) 
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Proof. As in the proof of Theorein 6.1 we can safely assume that ^ is a tree T 
with root r €V . We extend any / : i7+ R to a function / on fi by setting 
f{r}y = 1 Vy) = f{r}y = 1 Vy ^ r, ry^ = 0). Using Theorem 6.1, we then write 

Var+(/) = Var+(/) < Var(/) 
< gap(T, r, East)-i 5^ M (cx,r Var,(/)) 

where the constraints {c^^T^xeT have been defined right after (26). 

Let us examine a generic term /i ^c^-.r Var2;(/)^ with x ^ r. Remember 

that Cx,T — Cx.T and moreover Cx^Tiff) = if ryj, = 1 for all y ^ x. Fur- 
thermore, for any r\ such that there exists y ^ x with % = 0, li^{r}x = 
1 1 {'nv}y¥=x) = P- In conclusion we have shown that 

M (cx,r Var, (/)) =M^+)/x+(c,,rVar+(/)) Va; ^ r (31) 

We now examine the dangerous term ^Cr-,r Varr(/)^ = ii {VBXr{f)^- Be- 
cause of the definition of / we can safely rewrite it as 

H (WarrU)) = M (X{3j/#r:r,^=0} Varr(/)) 

Let us order the vertices of the tree T starting from the furthermost ones 
by first assigning some arbitrary order to all vertices belonging to any given 
layer (= same distance; from the root) and then declaring a; < y iff either 
d{x, r) > d{y, r) or d{x, r) = d{y, r) and x comes before y in the order assigned 
to their layer. Next, for any 77 such that r/y — for some y r, define 
^ = mm{y : % = 0} and let 7^ := {2 € T : z > ^} (see Fig 3). 




Notice that the siibgraph 7^ is again a tree and we define its root to be 
the ancestor t) of ^ in T. Then, using convexity of the variance, we can write 
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M(x{3y#r:,,„=o} Varr(/)) = M (x«#rM (Varr(/) | C)) < A* (xc#r Varr^ (/)) 

In order to bound from above Varr^ (/) we apply the Poincare inequality in 
7{ with constraints {cz^r^} and root v together with Theorem 6.1: 

Varr,(/) < gap(Z, East)"! ^ iit, {cz,r, Var,(/)) 

Notice that, by construction, Cz^T^irj) = Cz,t{v) for any z & T^, including the 
root w of 7^ where Cv^T^{r]) = 1 by definition and Cy^r{ri) = 1 because 77^ = 0. 
Putting all together we conclude that 

M (X{3 3/7^r: 7,^=0} Varr(/)) < gap(Z, East)"^ ^ IJ'{X{^vjir:r,v=0}Cx,rVaX:o{f)) 

< gap(Z,East)-V(-^+) ^ (c,,t Var+(/)) 
where we have used once more the observation before (31) to write 

If we now combine (6), (31) and (6) together we get 

Var+(/) < 2 (gap(Z,East)/x(r2+))"' ^ /x+ (c,,r Var+(/)) 

and the proof is complete. □ 
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